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We propose an alternative formalism to simulate CMB temperature maps in ACDM universes with 
nontrivial spatial topologies. This formalism avoids the need to explicitly compute the eigenmodes 
of the Laplacian operator in the spatial sections. Instead, the covariance matrix of the coefficients 
of the spherical harmonic decomposition of the temperature anisotropics is expressed in terms of the 
elements of the covering group of the space. We obtain a decomposition of the correlation matrix 
that isolates the topological contribution to the CMB temperature anisotropics out of the simply 
connected contribution. A further decomposition of the topological signature of the correlation ma- 
trix for an arbitrary topology allows us to compute it in terms of correlation matrices corresponding 
to simpler topologies, for which closed quadrature formulae might be derived. We also use this de- 
composition to show that CMB temperature maps of (not too large) multiply connected universes 
must show "patterns of alignment" , and propose a method to look for these patterns, thus opening 
the door to the development of new methods for detecting the topology of our Universe even when 
the injectivity radius of space is slightly larger than the radius of the last scattering surface. We 
illustrate all these features with the simplest examples, those of flat homogeneous manifolds, i.e., 
tori, with special attention given to the cylinder, i.e., T 1 topology. 
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I. INTRODUCTION 

It is becoming widely recognized that universes 
with a nontrivial spatial topology may be more natu- 
ral models for our Universe than the traditional sim- 
ply connected ones. This naturalness can be invoked 
from the mathematical point of view by arguing 
that there is an infinity of locally homogeneous and 
isotropic multiply connected 3-spaces, while there 
are only three simply connected ones; or with phys- 
ical arguments coming from incursions into the no- 
body's territories of Quantum Gravity and Quantum 
Cosmology. On the other hand, from a more prag- 
matic point of view, we can argue that cosmological 
models with a nontrivial spatial topology offer a very 
rich field of research, and are particularly well suited 
to explain certain reported "anomalous" features in 
CMB temperature maps, such as the alignments of 
their low ^-modes 0, 0] ■ 

Conversely, the full sky CMB temperature maps 
produced by the space missions COBE and WMAP 
provide us with an amazingly rich and high qual- 
ity amount of data with which we can look for the 
topology of space. This is very compelling for those 
who wish to unmask our Universe and see its shape, 
since Cosmic Topology is at present an almost ex- 
clusively observational and phenomenological issue, 
due to the lack of an accepted fundamental physi- 
cal theory which can predict the global topology of 
space. 
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Theory demands topology of space to leave sev- 
eral different kinds of marks in CMB temperature 
maps. Two of them have been largely studied and 
exploited to try to unveil the shape of our Universe, 
the distorsion of the angular power spectrum with re- 
spect to that of a simply connected universe 0-0, 
and the existence of "circles in the sky" |j| . Two 
other closely related sig natures, a non-null bipolar 
power spectrum 0, lllj . and alignments of the low 
£-modes 0, 0, have been only marginally used. 
Our main motivation for deciding to adventure into 
cosmic topology with the CMB was the desire to get 
a deeper understanding of the nature and properties 
of these alignments as a topological signature. 

One indispensable tool for a project like this is a 
software facility to produce simulated CMB temper- 
ature maps in multiply connected universes, so that 
we could systematically study the effects of differ- 
ent sizes and topologies on these alignments. These 
simulation procedures exist and have been used in 
several studies in cosmic topology E3~E3i so 
we could expect that this issue of the project would 
not present any problem. 

However, almost all known methods for comput- 
ing CMB temperature anisotropies in multiply con- 
nected universes need to solve the Helmholtz equa- 
tion in the manifold modelling our 3-dimcnsional 
space, the only exception to our knowledge being 
the work of Bond et al. ^3 ■ To solve the Helmholtz 
equation is a relatively easy problem in Euclidean 
manifolds H, fl(| , but a very difficult task in spher- 
ical [l8l Il9|| and hyperbolic 3-spaces [2(j. Indeed, 
the spherical case has been completely solved analit- 
ically only recently by Lachieze-Rey while for 
compact hyperbolic manifolds the only possible ap- 
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preach is numerical 20] . 

Among other results, in this paper we present 
a new approach to the computation of the corre- 
lation matrix (ae m a^, m ,) of the coefficients of the 
spherical harmonic decomposition of CMB tempera- 
ture anisotropies in a universe with nontrivial spatial 
topology. The main feature of this approach is that 
it avoids the explicit computation of the solutions 
of the Hclmholtz equation in the spatial sections of 
spacetime. Instead, we express the correlation ma- 
trix in terms of the covering group alone. Inciden- 
tally, the idea of generating a CMB map exploiting 
the symmetries of the quotient space was already 
sug gest ed time ago by Janna Levin and collaborators 
[l2l Il3 | . In particular we wish to quote a citation of 
a nontrivial claim in |l3j (p. 2695) which we have, 
in our opinion, succeeded in achieving: "By under- 
standing the symmetries of the fundamental polyhe- 
dron and the identification rules, a CBR pattern can 
be deduced without the need to explicitly obtain the 
spectrum mode by mode." 

Our main formal result is a generic decomposition 
of the form 

x r = x s - c - + x Ls - , (1) 

where X may be any covariance function which can 
be related to the two-point correlation function of 
the Newtonian potential (see Sec. IIII|) . T is the cov- 
ering group of (the multiply connected) space, s.c. 
stands for "simply connected" , and t.s. means "topo- 
logical signature" . 

Thus, in Eq. 0, X r is the covariance function 
computed in the manifold M = M/T, and X s - C - the 
same covariance function but computed in the uni- 
versal covering space M. It means that all the topo- 
logical information is encoded in the "perturbative" 
term X t s -, and that is why we refer to it, generically, 
as the topological signature of X r . This decomposi- 
tion is always possible whenever one can express X T 
in terms of T, as for example in the PSH method for 
detecting multiple copies of standard candles |2l| . In 
the present case we succeeded in writing the corre- 
lation matrix of the a^ m 's in this way by formally 
manipulating the two-point correlation function of 
the Newtonian pot ential derived by Bond, Pogosyan 
and Souradeep \u\ . 

Our approach to compute the correlation matrix 
(»i m ii(, m ,) has some methodological advantages in 
the simulation of CMB temperature maps. Indeed, 
by means of a suitable decomposition of the cover- 
ing group r in cyclic subgroups, we are able to write 
down a formula for the correlation matrix of a com- 
plicated topology in terms of the correlation matri- 
ces of the cyclic topologies (topologies with a cyclic 
covering group) that cover it maximally. Since cor- 
relation matrices for cyclic manifolds are relatively 
easy to compute (we obtain a closed quadrature for- 
mula for the cylinder), we expect to obtain in the 
near future more efficient ways to simulate maps for 
complicated manifolds. 



The decomposition of T in cyclic subgroups de- 
scribes in a transparent way the symmetries of the 
manifold, and this fact gives rise to another advan- 
tage of our approach, in this case from the observa- 
tional side. Universes with a cyclic topology present 
an alignment along the "direction of the generator 
isometry" . It follows that a CMB map in a universe 
with a nontrivial topology might present "patterns of 
alignment" (one for each ^-mode) characteristic of its 
shape and size. Indeed, and this shall become "obvi- 
ous" , symmetries of the quotient manifolds translate 
into symmetries of their patterns of alignment. We 
propose a method to search for these patterns by con- 
structing maps of the dispersion of the squares \ag m \ 2 
around the power spectrum. This opens the door 
to the development of methods to look for topology 
by searching these patterns, instead of limiting our- 
selves to considerations concerning only the special 
directions defined by the alignments. 

For the sake of brevity, other advantages of our 
approach to compute the matrix {aimO>% m t) ar e dis- 
cussed in Sec. I VII only. We prefer now to make a few 
remarks on some limitations of our work. We have 
considered here a few simplifications to develop the 
formalism, and worked out the details for the very 
simplest nontrivial topologies. In fact, we (i) have 
considered the Sachs- Wolfe effect as the only source 
of temperature anisotropies in CMB maps, (ii) have 
written the correlation matrix {atmO-% m i) formally 
only for Euclidean 3-spaces, (iii) have worked out 
the details for homogeneous flat manifolds, (iv) have 
made a detailed analysis and some simulations only 
for cylinders (T 1 topology), and (v) these simula- 
tions were done using the Einstein-de Sitter model. 
We wish to close this introduction by justifying each 
one of these simplifications. 

The shape of space is a global property, thus we 
expect the topological signatures in CMB to show 
themselves on very large scales only. Although a 
proof is missing, we believe that the main features of 
these signals would observationally appear if we re- 
strict the searches to the low £-modes in the temper- 
ature maps, i.e., we do not need high resolution CMB 
maps in Cosmic Topology! Since the Sachs- Wolfe ef- 
fect is the main source of temperature anisotropies 
at these scales 01 > we expect that theoretical ex- 
plorations considering only this effect will put in ev- 
idence the main features of the topological signals 
that we would observe in a real map. The addition 
of the missing part of the anisotropies will only mod- 
ify quantitatively the predictions made with our ap- 
proximation, and thus, will only be important when 
adjusting theoretical models with data. 

We consider this paper technically hard, so much 
care has been taken to write it in a clear and ped- 
agogical way. The main features of our formalism 
and of the topological signatures we predict in CMB 
can be understood by restricting the presentation to 
Euclidean topologies. The inclusion of nonzero spa- 
tial curvature will only introduce additional technical 
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considerations (and nothing qualitatively new), thus 
we decided to leave the spherical and hyperbolic cos- 
mological models for a future paper. 

The same applies to the lack in the paper of de- 
tailed explicit correlation matrices for nonhomogc- 
neous flat manifolds. Indeed, the price we pay for 
simplicity and transparency in the presentation of 
the results for each specific cyclic topology is the 
need for very hard calculations in the middle steps, 
as can be seen in the appendices. Explicit calcu- 
lations for nonhomogeneous flat manifolds will only 
add one page to the main body of the paper, and one 
or two more appendices to the already large list of 
them (see the end of the introduction). An exhaus- 
tive presentation for all the Euclidean manifolds is 
left for a future work. 

The last but not the least, we performed the sim- 
ulations with the Einstein-de Sitter model for sim- 
plicity. Nothing is lost from the theoretical point 
of view with this simplification since, as discussed 
in the paper, the structure of the topological signa- 
tures in CMB is captured in this oversimplified and 
old fashioned model of our Universe. However, more 
realistic ACDM models will be required to confront 
theory with observations quantitatively. 

We close this introduction by giving a detailed de- 
scription of the structure of the paper. In Sec. ITTlwc 
briefly review the two most common methods to sim- 
ulate CMB temperature maps in multiply connected 
universes, as well as present the method we have de- 
veloped, and implemented for the Euclidean case. In 
Sec. IIIII we define the topological signature in a cor- 
relation function, perform the decomposition of the 
covering group of a quotient space in its cyclic sub- 
groups, and write the topological signature in terms 
of this decomposition. We also show here that the 
symmetries of a quotient space appear transparently 
in the decomposition of its covering group in cyclic 
subgroups. 

In Sec. II VI we apply our formalism to the homoge- 
neous Euclidean manifolds, which are the simplest. 
We first compute the correlation matrix and the an- 
gular power spectrum for the cylinder, and apply 
the general results of the previous section in order 
to write down the correlation matrix of the spherical 
harmonic coefficients and the angular power spec- 
trum for a generic torus. Finally, as an illustration, 
we write those expressions explicitly for the chimney 
(T 2 topology). 

In Sec. |V|we first show, by means of simulations, 
that in a universe with the topology of a cylinder the 
low ^-modes are aligned in a similar fashion as they 
are in the WMAP data. We then use the results in 
the previous sections to argue that CMB tempera- 
ture maps in a universe with a nontrivial topology 
must present characteristic patterns of alignment, 
and propose the method of mapping on the sphere 
the dispersion of the squares |<z^ m | 2 to look for them. 
Finally, in Sec. I VII we discuss in detail the results of 
this paper and suggest further lines of research. 



In brief, the main goal of this paper is to show 
that our approach to the computation of the corre- 
lation matrix (ag m ag, m ,) (Sec. [HJl, together with the 
decomposition of the covering group of a manifold 
in cyclic subgroups (Sec. IHIfl . led to the discovery 
of a new topological signature in CMB temperature 
maps, i.e., the "patterns of alignment" (Sec. EJ. To 
illustrate this we have used the simplest example, 
i.e., that of fiat homogeneous manifolds (Sec. IIV(I . 

The paper has four appendices. In Appendix lAl 
we collect standard definitions and results related 
to spherical harmonics in order to set the conven- 
tions used in this paper. The technical calculations 
needed for the computation of the correlation matrix 
(a^ m a^, m ,) for the cylinder are presented in Appen- 
dices[B]and[U] In the former we develop from scratch 
the theory of Clausen ^-functions, for which we have 
not found any suitable reference in the literature. In 
the latter we compute a function that is the key part 
for computing efficiently the topological signature of 
the correlation matrix for the cylinder. Finally, in 
Appendix iDl we briefly reproduce known results for 
the correlation matrix of five out of the six compact 
orientable Euclidean 3-spaces. These formulae have 
been obtained previously in the literature by consid- 
ering explicitly the solutions of the Hclmholtz equa- 
tion in these manifolds, and are written in terms of 
the A;-modes jj, Il6| . Our derivation avoids the need 
for considering these solutions. 

II. SIMULATING CMB TEMPERATURE 
MAPS 

In this section we briefly describe two methods 
currently available for simulating CMB temperature 
maps in universes with non-trivial spatial topology, 
and proceed to develop our own formulation. We 
consider ACDM universes, where the background 
metric of spacetime is of the Robertson- Walker type, 
and include scalar and adiabatic perturbations as the 
seeds for the temperature anisotropies of the CMB. 
In the Newtonian gauge we have 

ds 2 = a 2 (ri) [(1 + 2$)rf?7 2 - (1 - 2<P)j lj dx l dx j ] 

for the metric, where r\ is the conformal time, a(r\) is 
the scale factor, $ is the Newtonian potential, and 

la ■= ( 1 + f {x 2 + v 2 + z 2 )^J &ij 

is the metric of the spatial section of the background 
with sectional curvature K = 0, ±1. 

The matter content consists of radiation (£l r ), 
baryonic and cold dark matter (fl m = ilb + fl c d m ), 
and dark energy in the form of a cosmological con- 
stant (^a)- Since we are interested on fluctuations 
on large angular scales, we make the assumption of 
instantaneous recombination and do not consider fi- 
nite thickness effects. The main contribution to the 



4 



temperature anisotropy observed at the direction n 
comes from the complete (ordinary plus integrated) 
Sachs- Wolfe effect 



ST, , 1 _ 



drj 



(2) 



d<S> 

drj 



where the index LSS stands for 'last scattering sur- 
face', the index for present time, and R(rj) is the 
comoving distance at instant rj between a photon, 
scattered at rj^ss, and the observer. 
The Newtonian potential is written as 



$(r/,x) = J dqF(n,q)t(q,x) 



(3) 



The temporal part satisfies the equation 

F"(r ] )+m(l + c 2 s )F / (r 1 ) + [2H'+ (4) 
(1 + 3c 2 s ) (H 2 -K) + c 2 s q 2 ] F( V ) = , 

where c s is the speed of sound in the fluid and Ti = 
a! I a is the Hubble parameter in conformal time. On 
the other hand, the spatial part consists of solutions 
of the Helmholtz equation 



(A + g 2 K(g,x) = 



(5) 



where the index q has been put as a variable in £ for 
simplicity of notation. 

The integral in Eq. has to be understood in 
a measure theoretic sense. Indeed, for multiply con- 
nected spaces the measure dq is not the usual one 
but a combination of a discrete and an absolutely 
continuous measures, reducing the integral in (pj to 
a sum and an integral in the usual sense. In partic- 
ular, if the space is compact, the measure reduces to 
a discrete one. This comes from the well-known fact 
that not every eigenmode of the Laplacian operator 
in the universal covering space M is also an eigen- 
mode in a quotient space M — M/T. In fact, only 
eigenmodes in M satisfying the invariance conditions 



(6) 



for any g € T project to eigenmodes in M. 

The most straightforward way of simulating CMB 
temperature maps is by solving Q and (J5J , perform- 
ing the sum in © , and then evaluating the SW effect 
(J2J. However, one has to consider that a temperature 
anisotropy map is a realization of a random field on 
the 2-sphere, and this randomness is inherited from 
that of the Newtonian potential ©. There are cur- 
rently two ways to implement this random character 
in the simulations, and one goal of this paper is to 
propose a third one. 

The first and most direct method is to consider the 
randomness in the temporal part of the decomposi- 
tion © of the Newtonian potential. The two-point 



correlation function of the Newtonian potential at 
fixed time 77 can then be written as 



($(7 7 ,x)<D(7 7 ,x')> = J dqdq'f( V ,q,q')x (7) 

?(g,x)fV,x'), 

where f(rj, q, q') — (F(r), q) F(r), q 1 )) is the two-point 
correlation function for the amplitudes of the scalar 
perturbation modes and £(g, x) are normalized solu- 
tions of the Helmholtz equation. 

Assuming that the Newtonian potential is a ho- 
mogeneous and isotropic random field, the two point 
correlation function (JJJ reduces to a function of 
time rj and the distance d(x, x'), and thus we get 
f(V) 9> l') = p <s>(v, q)S(q - q'), where P^>(r], q) is the 
gravitational power spectrum. If, in addition, the 
Newtonian potential is assumed to be gaussian, its 
random character is completely encoded in the vari- 
ance of the temporal part 



(F 2 (ri, q)) = P<s>(r], q) 



(8) 



Specifying this function, one then takes as an ini- 
tial condition, F(r) init , q), a realization of a normal 
distribution with zero mean and variance given by 
Eq. (jSJ), and some suitable condition for the ini- 
tial first derivative. With these initial conditions one 
solves for (@J , so one can now compute the potential 
(0. Topology is considered by restricting in © to 
normalized solutions of the Helmholtz equation sat- 
isfying the invariance conditions ©• This method 
has been extensively used in 0, 0, although 
in the former the authors do not consider the ran- 
domness of the function -F(?7, q). Instead the random 
character of the CMB maps is attributed exclusively 
to the random character of the eigenmodes of the 
Laplacian in compact hyperbolic spaces. 

The second method to produce simulated maps 
of CMB temperature anisotropics in universes with 
nontrivial spatial topology was first described in 0] , 
and used in [THl Ilr3 |. It is based in considering 
the randomness of the Newtonian potential in the 
eigenmodes of the Laplacian operator. We begin by 
decomposing the general solution of the Helmholtz 
equation, in the universal covering space, as a sum 
of fundamental solutions 



£(<?: x ) = ^2&rn(q)yim(q,X) 



£,m 



where 



3^ m (g,x) = p£ (q,x)Y im (n) 



(9) 



(10) 



is the normalized solution of the Helmholtz equation 
after separating it in radial and angular variables. 
Here we have put x = |x|, n is the unit vector in 
the direction of x, and the Y^ m (n) are the spher- 
ical harmonic functions (see appendix 0. Since 
the solutions l|l(J|) are normalized, the randomness of 
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the eigenmodes' amplitudes relies on the coefficients 

Introducing © in J3J, and using © and lfTU|) . 

we arrive at the decomposition of the temperature 
anisotropy map in spherical harmonics 



ST 



with multipole coefficients 



dqli m (q)Gt(q) , 



and the effects of physical cosmology given by 



(11) 



(12) 



Gi{q) = ^F(T] LSS ,q) pe(q,RLss) 



rrio dF 
2 drj 

JriLSS 



drj 



(13) 
pe(q,R(r})) . 



At this point it is convenient to recall how does 
topology enter in the story. Note that, due to the 
invariance conditions (jSJ, not every solution of the 
form jnj) is a solution in a quotient space. However, 
one would expect that any solution in a quotient 
space could be written in this form, the only problem 
being to find the correct coefficients £,e m (q)- These 
coefficients are not independent one from the other, 
since the invariance conditions © establish certain 
relations among them. In what follows we will as- 
sume that these relations can always be found, so 
that we will always represent an eigenmode in a quo- 
tient space by Eq. Q, with suitable coefficients. 

A crucial point in [lj] is the decomposition of these 
coefficients as 



itm{q) = \J P$(q)e£m(q) 



(14) 



where P®(q) is the gravitational initial power spec- 
trum, and the eg m (q) form a multivariate gaussian 
random variable, with a non-diagonal covariance 
matrix due to the relations among the coefficients 
£,Zm{q) coming from the invariance conditions. 

The simulation procedure can now be described. 
First we solve Eq. Q using the initial condition 
F{rji n it, q) = 1, and a suitable condition for the first 
derivative, and use this in (|13|) to compute Ge(q). 
Then generate a realization of the random variable 
e£ m (g) and use (|14|l to compute £,t m {q)- The map 
is now simulated by computing the coefficients ag m 
using l|12|) . and performing the sum in (Q) . 

An alternative method of simulation, also pro- 
posed in ^4|, is to construct the covariance matrix 
of the a£ m 's as 



(a hn a* e , m ,) = J dqdq'Gt>{q)Gt>,{q') x (15) 

& m {q)€i> m >{q')) ■ 

The substitution of (|14fl into l|15|l . and the evalua- 
tion of the resulting integral give rise to expressions 



for the covariance matrix in terms of the eigenval- 
ues and eigenmodes of the Laplacian operator. The 
multipolar coefficients are then obtained directly as a 
realization of a gaussian distribution with zero mean 
and covariance given by (|15|) . 

The method we propose in this paper lies along 
these lines, but we are able to manipulate the cor- 
relation function for the {;i m (q) in a way that avoids 
the need for an explicit determination of the eigen- 
modes of the Laplacian. Instead, the final expression 
after the integration of (|15fl is given in terms of the 
isometries of the corresponding covering group. 

Our starting point is an expression, obtained by 
Bond et al. in [llfl, that relates the two-point cor- 
relation function of the Newtonian potential in a 
simply connected universe, and that in a multiply 
connected universe, when both potentials have the 
same initial power spectrum. For a homogeneous 
and isotropic random field the expression is 



{$( V , x)$(r?,x')) r = \g\ ($(77,x)$(r,,. 9 x')) s - c - , 

(16) 

where \g\ = 1 if g is orientation preserving, and — 1 
otherwise. 

We now show how to use Eq. I|16(l in order to ex- 
press (II 51) in terms of the covering group, and for 
simplicity we will restrict the presentation to flat 
topologies. In Euclidean space, the most general so- 
lution of Eq. Q is written in the form 



£(g,x)= / d 3 kS(q-k)^)< 



(17) 



If we now expand the plane waves in spherical har- 
monics as 

e lkx = i e j t (kx)Y; m (n ]l )Y tm (n) , (18) 

£,m 

where je(x) is the spherical Bessel function of order 
£, and introduce it in l|17|l we obtain x) expressed 
as in Eq. © with 

l, m {q) = 47^ j d 3 k S(q - k) f(k)y/ m (n k ) , (19) 

where n k is the unit vector in the direction of k, and 
Pt{q,x) = jt(qx). 

Note that we have not decomposed £e m (q) as m 
Eq. JHJ. Instead, the decomposition i|19|) allows 
us to implement the randomness of the Newtonian 
potential in the modes £(k). In fact, introducing 
<|19fl in (|15fl . the covariance matrix for the aim's now 
reads 

ifl tm a% m ,) = (4ir) 2 t e - e '[d 3 kd 3 k'G e (k)G i ,{k')x 



(C(k)r(k0)y/ m (n k )^ m '(n k 



(20) 



G 



It is the correlation function of the modes £ (k) that 
carries all the topological information, as we will see 
in the following. 

Introducing l|17|) in J3J we obtain 

$(tj,x) = J d 3 kF(7],k)^(k) e ik x , 

thus the two-point correlation function of the New- 
tonian potential now reads 

($(r?,x)$(r?,x')) = J d 3 kd 3 k'F(r 1 ,k)F(r 1 ,k r ) x 

(£(k)£*(k'))e l ( k - x - k '- x ') .(21) 

At this point we have to recall that an Euclidean 
isometry can always be written as g = (R, r) , where 
R is an orthogonal transformation and r is an Eu- 
clidean vector, and that this isometry acts on a vec- 
tor x as gx = Rx + r. It is now an easy task to 
deduce from Eqs. JTEJ and J5lJ that 

<e(k)r(k')) r - ^(c(k)r(i?k')) s - c e - iflk '- r . 

In most inflationary models the initial perturba- 
tions of the gravitational field are homogeneous and 
isotropic Gaussian random fields, thus the correla- 
tion matrix of the k-modes in a simply connected 
universe takes the form 

(c(k)r(k')) s - c = ^«5(k-k'). 

The use of <|2"2"|) now yields 

(f(k)f*(k')) r = ^ ]>>(k-i?k') e -^ k '- , 

for the correlation matrix of the k-modes in the quo- 
tient space M = M/T, which when substituted in 
l(2"n)l finally gives 

(a lm a* e , m ,) r = (4^) 2 i 1 - 1 J — *„,(*) x (23) 

T^ m /(k)y^(n k ) , 

where the physical effects are encoded in 

\Mfc)=P*(fc)G*(fc)G*(fc), ( 24 ) 
and the topological information in 

TLW=^e- ik ' r y <m (n rk ). (25) 
ser 

The integration in 12/il) is over the whole k-space. 
The topological information is carried in Eq. I|25|l . 
which automatically selects the eigenvalues of the 
Laplacian operator in M, This can be seen in Ap- 
pendix [D] where T^ m (k) is expressed in terms of 
Dirac's delta functions centered in the eigenvalues of 
the Laplacian operator in the corresponding quotient 
spaces. 



III. DECOMPOSITION OF V IN CYCLIC 
SUBGROUPS 

In this section we develop some formal results in 
order to proceed further. Especifically, we define the 
topological signature of any covariance function that 
can be decomposed as 

as for example, the two-point correlation function of 
the Newtonian potential and the correlation matrix 
{ai m a £' m '}- Then we work out a suitable decompo- 
sition of a covering group in their cyclic subgroups, 
and write down the topological signature in terms 
of this decomposition. We also show here that the 
symmetries of a quotient space appear transparently 
in the decomposition of its covering group in cyclic 
subgroups. It follows that the main result of this 
section is the ellucidation of how these symmetries 
manifest themselves in the topological signature of 
CMB temperature anisotropy maps. 

Let us begin by writing the obvious decomposition 

X r = X sc +X r , (26) 

where T = T \ {id}. The second term in the right 
hand side is the topological signature in the covari- 
ance function. The expressions we present in the 
following are analogous to l|26|l and are also rather 
obvious. 

It is convenient to introduce a notation, so natu- 
ral, that has been used in l|26|) without any previous 
definition. Let S be any subset of isometries of the 
covering space M, then a superscript S in the covari- 
ance function means 

x s = J2 x9 ■ 

ges 

Then if M = M/T is a quotient space and Ti C T 
is any subset of the covering group, we can immedi- 
ately write X T = X Vl + A r \ Fl . This expression is 
nothing but the simplest generalization of Eq. (|26() , 
which corresponds to the trivial case Ti = {id}. We 
get a further generalization as follows, let Ti and T2 
be any two subsets of the covering group T, such that 

ri n r 2 - r 3 , then 

X r = X Fl +X T2 - X Ti + X r \ (r iur 2 ) _ ( 27 ) 

We can now write the formal result we are interested 
in. Consider the subsets G\,...,G n C V such that 
for any i ^ j, Gi n Gj = H, then by induction on 
$Z7\ we get 

n 

A r = ^A G *-(n-l)A H + A r \ G , (28) 

i=l 

where G = UGj. 
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To move forward, let G\ = (gi) and G2 = (32) be 
two cyclic subgroups of B, and let 6 M be a lift 
to M of the position of the observer in M . We will 
say that gi and g 2 are conjugate by an isometry that 
"does not move the observer" if there exists an isom- 
etry r\ fixing and such that g\ — r\~ x g2"(). Note that, 
as a consequence, we have that d(0, giO) = d(0, g 2 0), 
where d(x,y) is the distance between two points x 
and y in M. By extension, we will also say that the 
groups Gi and G 2 are conjugate by an isometry that 
does not move the observer. In addition, we will say 
that g\ is a minimal distance generator of G\ (with 
respect to the observer) if d(0,giO) < <i(0,70) for 
any other generator 7 G G\. 

Now consider the isometries gi , . . . , g n € T that 
generate the cyclic groups Gi = (gi), and such that 
iiij^j then G t n G 3 = {id}. By using Eq. we 
immediately obtain that the topological signature of 
the covariance function can be decomposed as 

n 

i=l 

In the following we will be particularly interested in 
the case where the g^s are minimal distance gener- 
ators of the Gj's, and the latter form a complete set 
of groups mutually conjugate by isometries that do 
not move the observer. 

Decomposing the topological signature further 
along these lines, let g±, . . . , g n , hi, ■ . ■ , h m £ I be 
minimal distance generators of the groups Gi = (gi) 
and Hj = (hj), and let G = UG; and H = UHj. 
Moreover, suppose that the G^'s and the Hj's form 
two complete sets of groups mutually conjugate by 
isometrics that do not move the observer, and such 
that G n H = {id} and d(0,giO) < d(0, hiO). Then 
the topological signature can be decomposed as 

n m 

x f = J2 + ]T X s - + x r \( Guff ) . 

i=l i=l 

We can proceed along these lines again and again, 
and obtain the following decomposition, in cyclic 
subgroups, of the covering group B, 

00 ki 

r = U U r « > ( 29 ) 

i=U=l 

where gij € T is a minimal distance generator of the 
cyclic group r^ , and such that 

1. For each i £ N, the set {IV, . . . ,1^} is a 
complete set of groups mutually conjugate by 
isometries that do not move the observer. 

2. If i ^ i', the sets U^ =1 Tij and U.^IVj have 
the identity as the only common element. 

3. If i < %', then d(0,gnO) < d(Q,g ¥ iQ). 



Then the topological signature of the covariance 
function can be written as 

1=1 j=i 

Thus, to compute the topological signature of the co- 
variance function for any multiply connected mani- 
fold, it is enough to know how to compute it for 
manifolds whose covering groups are cyclic groups. 
It is now obvious that this decomposition will be par- 
ticularly useful for calculating the correlation matrix 
of the dim's for any compact manifold once we know 
how to calculate it for cyclic flat (twisted cylinders), 
spherical (lens spaces), and hyperbolic manifolds. 

Let us now show that the decomposition (|29|l de- 
scribes transparently the symmetries of the quotient 
manifold M — M/T. Actually, this decomposition 
contains the symmetries of the Dirichlct fundamen- 
tal polyhedron of M centered at the observer's po- 
sition € M, which is what one expects to recon- 
struct with cosmological observations. Recall that 
the Dirichlet fundamental polyhedron centered at 
e M is the set V C M defined by (see H2) 

V = {x e M : d(0, x) < d{gO, x) for any g € T} . 

The first thing to be noted is that, although the 
whole covering group enters in this definition, it 
turns out that, in order to effectively construct the 
Dirichlct polyhedron, we only need the minimal dis- 
tance generators (and maybe the first few positive 
powers) of the first few cyclic groups By and their 
inverses. 

In fact, for each g g T consider the semi-space 

H g = {xeM : d(0,x) < d(gO,x)} . 

Then it is obvious that the Dirichlct polyhedron is 
the intersection of all of these semi-spaces. How- 
ever, there is a high redundancy here, since for a 
sufficiently large positive power n, we may have 

n-l 

n c > 

k=l 

and so this and further powers of g^ do not effec- 
tively contribute to the polyhedron T>q. Addition- 
ally, if some H g effectively contributes to the poly- 
hedron, so does H g -i, thus the same argument holds 
for the inverses of the minimal distance generators. 
Moreover, note that due to condition 3. above, for 
a sufficiently large i, it may be the case that the 
semi-spaces H 9i . do not contribute effectively to the 
polyhedron. 

The faces of the Dirichlct polyhedron are subsets 
of the boundary planes of the semi-spaces effectively 
contributing to it. In fact, for each H g effectively 
contributing, the corresponding face is orthogonal to 
the geodesic joining and gO, and cuts it at its mid- 
dle point. It follows that the decomposition l|29|l de- 
scribes the symmetries of the Dirichlet fundamental 
polyhedron of M centered at the observer. 
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IV. FLAT HOMOGENEOUS MANIFOLDS 

We have seen in the previous section that, to com- 
pute the topological signature of CMB temperature 
anisotropics in a given manifold, we just need to 
know how to compute it for the cyclic manifolds that 
cover it maximally. In the flat orientable case the 
cyclic manifolds are twisted cylinders, i.e., manifolds 
with covering group generated by a screw motion. 
We will now focus on the simplest case, the fiat ho- 
mogeneous manifolds, which are generated by trans- 
lations only. 

The flat homogeneous manifolds are 3-tori or T 3 
manifolds (generated by three linearly independent 
translations), chimneys or T 2 manifolds (generated 
by two linearly independent translations) , and cylin- 
ders or T 1 manifolds (generated by one transla- 
tion). Thus, we first compute the correlation ma- 
trix {o-imO*£' m r) T for cylinders, and then show how 
the decomposition ll.'SOIl is used to compute the topo- 
logical signature in this matrix for two- and three- 
dimensional tori. We also show that the computa- 
tion of the angular power spectrum in tori is greatly 
simplified by this decomposition. 

A. The cylinder 

Let us consider a cylinder orthogonal to the z- 
direction, that is with covering group generated by 



the translation g = (I, a), with a = Le z , where dis- 
tances are measured in units of the radius of the last 
scattering surface Rlss- This choice of the coordi- 
nate system is very convenient since here the cylin- 
der appears invariant under (i) arbitrary rotations 
around the z-axis, (ii) the parity transformation, 
and (iii) the reflection on the y — plane. Thus, 
according to Sec. IA 31 we will end with a real corre- 
lation matrix with no m-dependent correlations and 
the multipoles t and I' correlated only when both 
are even or odd. 

It is convenient to recall here that our cylinder has 
injectivity radius equal to L/2, thus cylinders with 
L < 2 are "small" and might present topological 
copies of discrete sources and/or circles in the sky. 
On the other hand, "large" cylinders, i.e., those with 
L > 2, have undetectable topology with the methods 
currently available |23| . 

The covering group of the cylinder is labeled by 
the integers as g n — (I,na), with n e Z. We then 
have from I125H that all the topological information 
is encoded in 

TL(k) = ^e-™ w y lm (n k ), 

and thus the correlation matrix of the ag m 's for the 
cylinder is simply 



(47T 



,2 -l-t 



(31) 



r 



To reduce this integral we may use any of the fol- 
lowing two identities, either 



~ink z L 



2ir^25(k z L- 2tt P ) , (32) 



pel 



or 



J2 e -ink,L _ 1 + 2 cos(nk z L) . (33) 

n£Z n=l 

The first identity is obvious since the left hand side 
is the Fourier expansion of the right hand side. This 
option yields a formula of the kind obtained in Ap- 
pendix [DJ which expresses the correlation matrix in 
terms of the eigenvalues of the Laplacian operator on 
the cylinder. The second one still uses a parametriza- 
tion in terms of the covering group, and thus can be 
used to isolate the topological signature. 

In fact, using Q33Jl to evaluate (|31|l . and integrating 
in spherical coordinates, we get 



(aim a £' m ') r — { a lrn 



)*■*■ + (a em a2, m ,) F , (34) 



where the simply connected part is as usual 

(aim a t'm') S C ~ C^'SuiSmm' , (35) 

with the (simply connected) angular power spectrum 
given by 

CT ■ = (4n) 2 / — 9 u (x) , (36) 
Jo x 

and the topological signature for the correlation ma- 
trix is given by 



(a em a* e , ml f = {^? S™, d(2) 5 mm > x (37) 

r — 

Jo x 



with 



F%(x) =2^2 dy C os(nxy)Vr(y)VF(y) , (38) 

n=l J -l 
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FIG. 1: Shapes of the integrands inside integrals of the type in 1)37(1 for I — £' — 2, and compactification scales L = 1 
and L — 2 in units of the radius of the last scatering surface. 



where V™{x) is the normalized associated Legendre 
function (see Appendix lA"|l . As expected, we have 
ended up with a real correlation matrix with factors 

A mod(2) and r) 

°W dim O mm ' . 

After evaluating the series in l|38|) , it turns out that 
Fp%,(x) is a piecewise continuous function. In fact, 
in each interval [27rg, 2w(q + 1)], it is a polynomial of 
degree (£ + £' + 1) in n/x. Indeed, the final result is 

F% (x) = £ Fu> 9) Q(x - 2*q) x (39) 

Q(2n(q + 1) - x) , 

where 0(x) is the Heaviside step function, and the 
form of the polynomial T^, (x, q) in the g-th interval 
of length 2ir is 

i+i' 

F%,(x,q) = 4^(-l) fe ^(0)x (40) 

fc=0 

g2k+i(q) y-J - <W ■ 

Here gk{l) are polynomials of degree k in q, and 
V^) m (0) is the k-th derivative of the polynomial 

V Wm (x)=Pr(x)PP(x) (41) 

evaluated at the origin. In Appendix [5] we present 
recurrence relations for the polynomials gk{q), and 
all the technical steps that take us from l(3"5|) to (|3~9"|) 
can be found in Appendix IU1 

The integrals appearing in the topological signa- 
ture l|37|l can be easily evaluated since Eqs. I|39|) and 
(|40|l allow an exact and very fast computation of the 
function F^,{x), and the integrands decay very fast, 
as illustrated in figs^andl^l hi this figures we have 
adopted, for simplicity, a scale invariant Einstein-de 
Sitter model, thus 

Vw(x) ocji(x)ji>(x) . (42) 

The nice behavior of the integrands in (|37|l is not a 
consequence of this particular choice of *£m'(x). Ac- 
tually, the integrand in (|57J) always decays very fast 



because (x) and (x) are both decaying func- 
tions, thus the evaluation of the topological signature 
for the cylinder is always very efficient. 

The computation of the topological signature of 
the power spectrum reduces to a simple integral. In 
fact we obtain 

cf = (47r) 2 ^^*«(|) 

with 

m—~l 

Using (|38|) to perform this sum, the Addition Theo- 
rem for Spherical Harmonics (see Appendix^} yields 
immediately 

$ = 2^? £^*u(j)<px{x), (43) 

where <pi{x) is the first Clausen (^-function given in 
Appendix [5] 

In figlHK we show the low £-modes of the topo- 
logical signature of the angular power spectrum of 
a cylinder, normalized w.r.t. C| c ', as a function of 
its size L. We can see that the topological signa- 
ture is typically much smaller than the cosmic vari- 
ance, even for small cylinders which have already 
been discarded observationally as candidates for the 
shape of our Universe because of the lack of antipo- 
dal matched circles in WMAP data 0, ||| ■ Thus it is 
apparent that the angular power spectrum is not a 
good indicator to look for topology in this case. 

The correlation matrix given by l|34 fl -l|4U |) corre- 
sponds to a cylinder for which the direction of com- 
pactification is parallel to the z-axis. The correla- 
tion matrix corresponding to a cylinder with a dif- 
ferent orientation can be easily obtained from the 
previous one by simply rotating the celestial sphere. 
Thus, parametrizing the rotations with Euler angles, 
if R(a, /3,7) € SO (3) takes the z-axis to the direc- 
tion of compactification of the cylinder, the topologi- 
cal signature of the corresponding correlation matrix 
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FIG. 2: Shapes of the integrands inside integrals of the type in l)^7|l for I ■■ 
and L — 2 in units of the radius of the last scatering surface. 



5, and compactification scales L = 1 



can be computed using the expressions ljA2(l - i|A4|l of 
Appendix ^ yielding 

/~ ^* \T Aim 1 — m)a N . ( a a\ 

\aima tm ,) R = e x (44) 

5^4m.i(/9)«C'm 1 (/5) X 
mi 

(a^mi a £' mi ) r j 

since (ae m a>*ti m i) sc ' is rotationally invariant. More- 
over, the 7 angle does not appear in this expression 
since R z {l) in (|A3|I does not move the z-axis, and 
(ai m a* f , m ,) T is invariant under such rotations. 

It should be noted here that, no matter its orien- 
tation, the cylinder is always invariant under parity, 
thus its correlation matrix will always conserve the 
factor 5™ e ° d{2) . On the other hand, the correlation 
matrix will remain real as far as we perform rota- 
tions with a — 0, since in this case we do not rotate 
the cylinder around the z-axis, and thus it remains 
invariant under reflections on the plane y — 0. How- 
ever, any rotation of the cylinder (other than one 
with p = 7t) makes it non-invariant under azimuthal 
rotations, thus the correlation matrix of an arbitrar- 
ily oriented cylinder has m-dependent correlations. 
All these features can be seen explicitly in 144|) . 



B. tori 

In order to calculate the correlation matrix of the 
0£ m 's for a two- or a three-torus we use the de- 
composition (|29[) of its covering group in cyclic sub- 
groups. Let Tij = (gij) be the covering group of 
the cylinder generated by the element € T, and 
let us write Li — d(0,gijO), gi = (J, Lje z ), and 
Ti = (g-i)- In the Euclidean case, the orientation 



preserving isometries that do not move the observer 
are rotations, thus let Rij £ SO(3) be the rotation 
taking e z to the unit vector along QijO. 

Using the decomposition 13U|) . we easily write the 
topological signature for the torus as a superposition 
of topological signatures of rotated cylinders. In fact, 

oo ki ^ 

{aim a* llm ,) T = J2J2( aem a *i'm>)% , ( 45 ) 

i=l 3 = 1 

where the correlation matrices of the rotated cylin- 
ders are written in terms of the Wigner £)-functions 
and Euler angles, according to (|44|) . as 

{«im«lm>)% = e i(m '- m)a " x (46) 

^ ( dm 1 (flij ) ^m'mi (Aj ) x 
mi 

/ * \f 4 

and (Pij , otij ) are the angular spherical coordinates 
of the vector gijO, and ki is the number of cylinders 
of size Li. Since any group of translations is invari- 
ant under parity, from Sec. IA 31 we know that the 
correlation matrix for a homogeneous flat manifold 
has always the factor <5^° d(2> , and this is evident from 
l|45|l . since it is just a sum of correlation matrices of 
cylinders. 

The power spectrum is rotationally invariant, thus 
from (|45|1 one can easily write down the expression 
for the topological signature of the power spectrum 
of the torus as a superposition of topological signa- 
tures of power spectra of cylinders, 

oo 

C? = £>C?'. (47) 
i=i 
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FIG. 3: Topological signature of the power spectrum of (a) a cylinder, and (b) a chimney with square base, for the 
first four ^-modes, normalized w.r.t. C s t ' c ' , and as a function of the scale of compactification L. Note that for each 
multipole one can have a suppression or an excess of power depending on the value of L. For typographical reasons, 
here we write F* instead of T. 



Let us consider a chimney with square base for 
the sake of illustration. It is convenient to orient the 
chimney so that its covering group consists of trans- 
lations in the horizontal plane. We take as generators 
of the covering group the translations g\ — (I, a) and 
g-2 = {I, b), with a = Le x and b = Le y . 

It is more convenient to reparametrize the cyclic 
decomposition as follows. Parametrize each cyclic 
subgroup by a pair of integer numbers (p, q) as G pq = 
(ff2Si)- Clearly, if the greatest common divisor of 
(p, q) is r, then 

G pq < Gesl , 

where '<' means 'subgroup of. Thus we must re- 
strict the labels to pairs (p, q) of coprime numbers. 

The only exceptions are when (i) p = ±1 and 
q = and viceversa, and (ii) when p = ±1 and 
q = ±1. Thus the first two complete sets of cyclic 
subgroups conjugate by a rotation are {Gi^Go^} 
and {Gi,!, G-^i}. In both cases the conjugation is 
performed by a rotation of 7r/2 around the z-axis. 
The compactification lengths of the corresponding 
cylinders are = L^.i = L and L\ t \ = = 
\/2L respectively. The Euler angles (/3, a) to rotate 
the corresponding cylinders from the z-axis to their 
orientation in the chimney, according to l|46|l . are 
(3 = 7r/2 in all cases, and ev^o = 0, ao,i = 7r/2, 
= 7r /4 and — 37r/4, respectively. 

To write the remaining complete sets of cyclic sub- 
groups conjugate by a rotation let us define, for a pair 



of coprime natural numbers (p, q), with p > q > 1, 
the groups 

G$ = G Pg = <<?!<??> , G$ = G^ p = ( 9 p 2 g^) , 
= G qp = (gWi) , G$ = G- P , q = (ghI P ) ■ 

The compactification lengths are all equal to L pq = 
y/p 2 + q 2 L, and the Euler angles (j3, a) to rotate 
the corresponding cylinders from the z-axis to their 
orientation in the chimney, according to (|46|l . are 
(3 = ir/2 in all cases, and 



^ pq 


= arctan — 
P 


^pq 


7T 

= 2 4 


^■pq 


a (2) 
pq 


77 (i) 

= — — 

2 pq 




= 7T — 





respectively. 

Let us denote by T pq the covering group of the 
cylinder with compactification scale L pq and oriented 
along the z-axis. Then, putting all this together, us- 
ing H45f) and (|46|l . and taking into account the invari- 
ance properties derived in Sec. IA 31 the topological 
signature of the chimney with square base is 

(at m a% m ,f = 6ZT £ d e mmi (tt/2) X (48) 

mi 

di, mi (ir/2)W%;-™ , 

with 



W%, mi = 2 ((a tmi a;, m /"> + (-l) m / 4 (a emi al m f^) + 4 £ cosma« (a tmi a% m f»« , (49) 

I 



where the sum in (p, q) is evaluated only for pairs of coprime natural numbers (p, q) such that p > q > 1 . 





The topological signature of the power spectrum 
of the chimney with square base is simply 

Cj = 2 (d 1 '° + cf M ) + 4 ^ C[ M . (50) 

(p.?) 

Since the topological signature of the power spec- 
trum of a cylinder converges quickly to zero as a 
function of the compactification scale (see figOt); it 
follows that the sum in l|5(J|l also converges quickly. 

Moreover, the topological signature of the power 
spectrum of the chimney is larger than that of the 
cylinder. This is so because the £-th mode of the 
topological signature of the angular power spectrum 
of the cylinder oscillates very slowly. Thus from (|50(l 
this signature is slightly higher in the chimney, as 
can be seen in figOb- Actually, this is a general result 
that holds for manifolds whose covering groups are 
not cyclic. 

V. PATTERNS OF ALIGNMENT 

The nondiagonal character of the topological sig- 
nature of the correlation matrix of the a£ m 's in multi- 
ply connected universes and their m-dependence are 
manifestations of their globally anisotropic nature. 
They manifest themselves in statistically anisotropic 
temperature maps, i.e., realizations of random tem- 
perature fluctuations for which mean values of func- 
tions of the temperature over ensembles of universes 



depend on the orientation [10|. 

In this section we analyze an expected conse- 
quence of the topology of space on the temperature 
anisotropies of the CMB that has not received the 
deserved attention up to the present, namely the ex- 
istence of preferred directions in space. We show 
that the decomposition of the topological signature 
of the correlation matrix of the a£ m 's in a universe 
with a complex topology in signatures corresponding 
to cyclic topologies, demands the existence of "pat- 
terns of alignments" along these directions. For the 
sake of simplicity we consider the Einstein-de Sitter 
model, thus from now on we will take l|42(l to perform 
all our calculations. 

We want to call attention to the existence of align- 
ments of the low £-modes of the CMB temperature 
maps in multiply connected universes. Indeed, in 
fig^]we show a low resolution temperature map sim- 
ulation for a cylinder with L — 2 (in units of Rlss)i 
together with the maps corresponding to the first 
four ^-modes. One can see that these ^-maps present 
alignments along the z-direction, which in this case 
is the unique direction of compactification of space. 

Similar alignments as those present in our sim- 
ulations have been reported as being observed in 
WMAP data, and have been attributed to a possible 
nontrivial topology of space with the shape of a cylin- 
der 0. These models have been quickly abandoned 
due to the lack of circles in the sky which should be 
present if the Universe were small P, . However 
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FIG. 5: Plots of the topological signatures of the squares (|a£ m j 2 ) r , normalized w.r.t. Cf, for low multipoles 
(2 < I < 5), in a universe with cylindrical topology aligned with the polar axis, and as a function of the scale of 
compactification L. For typographical reasons, here we write Y* instead of Y. 



our simulations show that even in universes slightly 
larger, and so not presenting such circles, these align- 
ments should still be observable. Thus whether these 
observed alignments are a consequence of a nontriv- 
ial shape of our Universe is still an open question 

We will show here that if our Universe had a 
nontrivial topology, its CMB temperature map will 
present characteristic patterns of alignment, even if 
its size is somewhat larger than the observable uni- 
verse. Moreover, from the observed patterns of align- 
ment, we might be able to reconstruct the shape of 
space. 

In fig|U we show the topological signature of the 
low ^-modes of the diagonal part of the correlation 
matrix of the c^ m 's, normalized w.r.t. Cf , for a 
cylinder oriented along the polar axis, as a function 
of the size L of compactification. It is apparent that, 
for a given £-mode, there are multipole coefficients 
for which their expected values are above the mean 
(the angular power spectrum), and others for which 
these expected values are below it. This is the reason 
why the low ^-modes in a cylinder are aligned. Ac- 
tually, the expectation values (|a^ m | 2 ) are all equal 
to Ci only in the simply connected case, thus the 
dispersion around the mean 



^ = w^ TT E((K™i 2 )-^) 2 ( 51 ) 

y m 

is null. However, in a multiply connected universe, 
this dispersion is non zero, and in a particular map, it 
adds to the cosmic variance. Thus it seems natural to 
propose the dispersion of the squares |a^ m | 2 around 



their mean value as a measure of these alignments in 
a map. 

In fig we show a plot of the dispersion l|51|) . nor- 
malized w.r.t. Cf, for a cylinder oriented along the 
polar axis, as a function of L, for low multipoles. 
Note that even for a large cylinder (L w 2) the dis- 
persion is larger than 15% of the power for multipoles 
up to £ = 5. Indeed, on these scales the dispersion is 
of the order of the cosmic variance, and thus might 
be detectable. 

In order to show that this is a good measure of the 
alignment of multipoles, and that it provides an effi- 
cient method to determine the directions of possible 
alignments in real or simulated maps, let us compute 




L 



FIG. 6: Dispersion of the (|a« Tn | 2 ) r , normalized w.r.t. 
C\ , for low multipoles (2 < I < 5), in a universe with 
cylindrical topology aligned with the polar axis, and as 
a function of the scale of compactification L. 
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FIG. 7: Dispersion of the (|fi^ m | 2 ) r , normalized w.r.t. Cf, for low multipoles (2 < I < 5), in a universe with 
cylindrical topology, as a function of its orientation relative to the polar axis, and for different values of the scale of 
compactification L. 



the dispersion of the squares (|a^ m | 2 ), Eq. J5TJ, for 
a cylinder which is oriented along a direction mak- 
ing an angle with the z-axis. Each one of these 
squares can be computed with 

(Ki 2 )« = E[C(«] 2 (Ki 2 ) r . 

mi 

which is nothing but l|44|) restricted to the diagonal 
part. 

In fig0it is shown this dispersion as a function of 
(3 for different multipole coefficients and for different 
values of L. One can see that the dispersion has a 
maximum when the cylinder is oriented along the z- 
axis. Thus in order to look for the alignments in a 
hypothetical universe with the shape of a cylinder, 
one should just rotate the celestial sphere around 
different directions until find those two opposite ones 
along which the dispersion of the squares \ag m \ 2 is 
maximum. 

However, in order to collect definitive evidence 
that the universe is indeed a cylinder, one should 
map the dispersion of the squares |a^ m | 2 on the 
sphere for each £-mode, i.e. one should determine 
the dispersion (|51J) as a function of the orientation 
of the celestial sphere. If the universe had the topol- 
ogy of a cylinder, these dispersion maps should be 
axially symmetric around a special direction, where 
the dispersion is maximum. Moreover, this direction 
should be identified with the direction of compacti- 
fication of the cylinder. 

If the universe has the topology of a flat homoge- 
neous manifold note, from l|45l) and (|46[) . that the 
topological signature is a superposition of rotated 



cylinders of different sizes. Thus a CMB map for a 
universe with this kind of topology might present 
alignments along the directions corresponding to 
these cylinders. In fact, rotating the celestial sphere 
and computing the dispersion of the squares |a^ m | 2 , 
an easy computation shows that whenever we per- 
form the rotation i?(0, —8, —</?), with 8 = flij and 
tp = otij, one has the cylinder labeled by (i, j) ori- 
ented along the polar axis, and thus dispersion maps 
might present local maxima along these directions. 

Whether these local maxima are observable in a 
given dispersion map will depend on (i) the scale 
of compactification of the corresponding cylinder Li , 
(ii) the background due to the simply connected part, 
and (iii) the other cylinders' topological signatures. 
For large values of Li the corresponding local max- 
ima will not be observable, however one can expect 
those maxima corresponding to the smaller cylinders 
to be detectable. The existence and distribution of 
these maxima in each dispersion map, together with 
their relative intensities is what we call a pattern of 
alignment. 

It might seem that the problem of constructing 
dispersion maps for manifolds that are not flat homo- 
geneous is more involved, since general cyclic man- 
ifolds do not have axial symmetry as the cylinder 
has. Eq. I|44|) depends on two angles only because 
the cylinder is axially symmetric, but in the general 
case the expression for the correlation matrix in a ro- 
tated frame depends on the three Euler angles. Thus 
it seems at first sight that in these cases, a dispersion 
map should be a function on the 3-sphere. Fortu- 
nately, the diagonal elements of the rotated correla- 
tion matrix depend only on the last two Euler angles 
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as 

{Wm\ 2 )l = E e i( "" l)7 x 

((3)d 

III III 2 

/ * \r 

thus the same conclusion holds in the general case. 
Dispersion maps on the 2-sphere for low £ -modes 
should display patterns of alignment showing the 
symmetries of our Universe if it has a (not too large) 
nontrivial topology. 

VI. DISCUSSION 

In order to study systematically the effects of a 
nontrivial spatial topology in the temperature fluc- 
tuations of the CMB, we need to have the ability 
to simulate efficiently temperature maps in multi- 
ply connected ACDM cosmologies. Almost all the 
usual methods to perform these simulations use ex- 
plicitly the solutions of the Helmholtz equation in 
3-manifolds with nontrivial topology. The compu- 
tation of the eigenfunctions and eigenvalues of the 
Laplacian operator is simple only in Euclidean man- 
ifolds, while in spherical and hyperbolic spaces it is 
a nontrivial problem. In fact, it is only recently that 
an analytical computation has been achieved for all 
the spherical manifolds. The hyperbolic cases still 
have to be done numerically. 

In this paper we have developed a simulation pro- 
cedure that avoids the explicit use of the solutions 
of the Helmholtz equation. Instead, our results are 
expressed in terms of the covering group T of the cor- 
responding manifold. In this section we summarize 
the details of the method, its efficiency, the simple 
applications performed here, and discuss future re- 
lated work. 



A. The formalism 

The cornerstone of our method is formula (TSJ, 
which is the two-point correlation function of the 
scalar perturbations in a multiply connected universe 
expressed in terms of the covering group of the man- 
ifold 0| . By means of simple formal manipulations 
we obtain an expression for the correlation matrix of 
the spherical harmonic expansion coefficients of the 
temperature maps, Eqs. I|23(l - H25|l . which contain all 
the topological information expressed as a sum over 
the elements of the covering group. 

Former applications of (|16f) required a regulariza- 
tion procedure in order to account for divergences 
of the series, as well as some resummation tech- 
niques for accelerating the convergence. We do not 
have these problems here because the divergent se- 
ries, which are actually distributions, appear only 
inside integrals. Indeed, on the one hand, we show 



in Appendix IdI that our formalism easily reproduces 
results previously reported in the literature, as well 
as some simple generalizations, without the need of 
any regularization procedure. On the other hand, 
elementary decompositions of the two-point corre- 
lation function 1)16(1 . shown in Sec. IIIII guarantee 
that our final expressions are highly convergent, as 
discussed below. 

Two decompositions of a generic covariance func- 
tion which can be written as a sum over the covering 
group are crucial for the efficiency of our formalism. 
The first one, a trivial decomposition given by (|2(j|) . 
defines the topological signature of the covariance 
function. When written for the correlation matrix 
of the harmonic expansion coefficients, it yields the 
topological signature in the temperature anisotropy 
maps, as illustrated for the cylinder by (|34|l ■ This ex- 
pression shows that the topological signature is noth- 
ing but a "perturbation" of the correlation matrix 
corresponding to the simply connected case. Since, 
as discussed in Sec. IIV Al these "perturbations" are 
small, the efficiency of the calculation follows. 

The second decomposition given by (|5U|l allows us 
to write the topological signature of any manifold 
in terms of the topological signatures of its maximal 
covering manifolds with cyclic covering groups. The 
example of the tori illustrates the power of this ap- 
proach, since we can write down explicit formulae 
for a general torus whether its generating transla- 
tions are orthogonal and/or equal. Trying to do this 
with the explicit use of eigenfunctions of the Lapla- 
cian (or with the method used in Appendix [DJ turns 
out to be tedious if not difficult. 

Moreover, by construction, this decomposition is 
invariant under the symmetries of the manifold, thus 
it carries information on how shall these symmetries 
manifest in individual CMB temperature anisotropy 
maps, as will be discussed in Sec. IVI Bl 

Another advantage of this second decomposition is 
the simplicity for writing down the power spectrum 
for complicated manifolds. Expressions like l|47l) and 
H5Ufl are computationally very efficient once we have 
saved the power spectrum for cyclic manifolds as 
a function of its scale of compactification, since we 
have just to perform a weighted sum of power spec- 
tra for cyclic manifolds at different scales considering 
the multiplicity of the decomposition. 



B. Topological signatures 

A further advantage of splitting the correlation 
matrix of the multipole coefficients into its simply 
connected part and its topological signature is that 
we can identify very easily the geometric features of 
the signature. Although we have made explicit cal- 
culations for flat homogeneous manifolds only, qual- 
itatively these results are general. 

Universes with cylindrical topology of size L w 2 
present clear alignments of their low ^-modes along 
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the direction of compactification. A dispersion map 
of the squares (|a£ m | 2 ), for a given low £, exhibits an 
axial symmetry around this direction, thus it reduces 
to a function of the polar angle. These dispersion 
maps are shown in figd 

By decomposing the covering group T in cyclic 
subgroups one can see that, whatever the shape of 
our Universe, and if it is not too large, dispersion 
maps (one for each individual low multipole) might 
show patterns of alignment. In the general case such 
maps are functions on the two-dimensional projec- 
tive space or, by a lifting, on the 2-sphere. Although 
we have shown the existence of patterns of alignment 
explicitly only for homogeneous flat manifolds, it fol- 
lows from the exposition of the general formalism 
that the same conclusions hold for any manifold of 
constant curvature. Thus, we propose the construc- 
tion of these dispersion maps in the WMAP data, 
and so the search for patterns of alignment, as a new 
method for detecting a possible nontrivial topology 
of our Universe. 

It is interesting to comment on some features re- 
lating Levin and collaborators' proposal of pattern 
formation in CMB temperature maps and the results 
we present in t his p aper. The patterns proposed by 
Levin et al. [T^-ll3j are due to individual eigenmodes 
(fc-modes), the patterns we have identified here are 
due to multipole modes (£-modes). In either case 
the modes compete to form their patterns in a CMB 
temperature map, however the observable modes in 
a map on the sphere are the latter, since spherical 
harmonics form a base on the space of functions on 
the sphere. On the other side, the association be- 
tween real space perturbations and angular temper- 
ature fluctuations requires some averaging over the 
fc-modes [(|, thus these patterns appear mixed in a 
map and their observation might demand more elab- 
orated techniques. 



C. Further remarks and future research 

The formalism we have developed in this paper 
reveals new insights on the problem of characteriz- 
ing the marks that topology leaves in CMB maps, 
and opens up new possibilities for developing further 
methods for unvealing the shape of our Universe. 
It makes explicit that the multipole alignments ob- 
served in COBE and WMAP full sky CMB tempera- 
ture maps may be a manifestation of its global shape, 
provides details of the nature and features of these 
alignments, and gives at least one methodology to 
test this hypothesis. As a consequence, further work 
is much needed. 

One line of further research is the implementa- 
tion of our formalism in the spherical and hyperbolic 
cases. One way to do this requires first to identify the 
radial part of the fundamental solution 1)11) [) of the 
Helmholtz equation in the universal covering, as well 
as the analog of the "plane wave expansion" solution 



<|17fl in these geometries, and to write the expansion 
of the corresponding "plane waves" in spherical har- 
monics as in <|18[) . The difficult part seems to be 
expressing the "plane wave expansion" in a suitable 
form to reproduce the formal steps used in the Eu- 
clidean case. 

Moreover we have to compute the topological sig- 
nature of all other cyclic manifolds, in order to ex- 
tend the computations to any quotient space that 
could be a candidate for the shape of space. We also 
need to include acoustic oscillations, and Dopplcr 
and finite width effects in ACDM models so that we 
could determine the relevant angular scales in Cos- 
mic Topology, i.e. the angular scales at which the 
topological signatures appear. This is a crucial step 
in order to confront quantitatively the theory with 
real CMB maps in an efficient and rigurous way. An 
ultimate goal may be to implement all this method- 
ology in a software package for public use. 

The identification of the "topological signature" 
of the correlation matrix (a£ m a*ii m i) T also opens up 
a path for solving a problem raised by Riazuelo et 
al. in 14] . The correlation matrix for a multiply 
connected universe is non diagonal and, tipically, in- 
dependent. In fact, this is the source of the statistical 
anisotropy of the CMB in these universes. However, 
for very large manifolds this correlation matrix be- 
comes effectively diagonal, and equal to that corre- 
sponding to the universal covering counterpart. A 
natural question raises up: at what typical scales 
does the correlation matrix "becomes diagonal" ? In 
terms of our formalism this problem can be stated as 
finding the scales at where the topological signature 
becomes observationally negligible compared to the 
simply connected part. A closed analysis of the topo- 
logical signature might give some answers to this and 
related questions. For example, establishing bounds 
on the integral in l|37fl might solve the problem for 
the cylinder. 
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APPENDIX A: SPHERICAL HARMONICS 

In order to be self contained and to set the nota- 
tion used in the paper, in this appendix we present 
basic definitions, some useful formulae of spherical 
harmonic functions and Wigner rotation matrices, 
as well as some invariance properties of the corre- 
lation matrix (a£ m a^, m ,) under coordinate transfor- 
mations. For a complete treatment of s phe rical har- 
monic functions the reader can consult 25] . 



A very useful formula is given by the Addition 
Theorem for Spherical Harmonics 

47T 

P,(n ■ n') = — - ]T F, m (n)Y7 m (n') , (Al) 

m=-t 

which for the particular case n = n' yields the iden- 
tity 



n2 2£+l 



1. Basic definitions 

Let us denote by n = (6, if) a point in a 2-sphere 
parametrized in the usual spherical coordinates, then 
the spherical harmonic functions are defined as 



Y tm (n) = J^±l ^ Pt(co S 9) e" n * 



47r (£ + m)\ 



where 



pp( x ) = (-ir(l- x T /2 £;Pt(x) 

are the associated Legendre functions with non- 
negative index < m < £, and with 



1 



1 



2 l £\ dx e 

being the Legendre polynomials. The associated 
Legendre functions with negative index m are de- 
fined by 

Moreover, it is often convenient to introduce the nor- 
malized associated Legendre functions 



'2£+l (£-m)\ 



(£ + m)\ 



PPM 



It can easily be seen that the Legendre polyno- 
mial Pt{x) is an t-th. degree polynomial of parity £, 
and thus the associated Legendre function PP(x) is 
a function of parity £ — m. It follows that the func- 
tion Tu' m {x) defined in (|TT1) is an (£ + £')-th degree 
polynomial of parity £ + £' , and thus the expression 
for F^,(x) in which is evaluated in Appendix 

lUl contains only even polynomials. 

Spherical harmonics form a complete orthonormal 
set of functions on the sphere, thus their most com- 
mon application is in the expansion of functions, like 
a CMB temperature anisotropy map, in multipoles 
as in Qllfl. where the coefficients ai m , called the mul- 
tipole coefficients, are given by 

dO^(n)F/ m (n) . 

Since the temperature map is a real function on the 
sphere, the multipole coefficients obey the constraint 



2. Wigner D— functions 

In several ocassions it is convenient to rotate the 
sphere and compute the multipole coefficients in this 
new coordinate system. This can be achieved by 
means of the Wigner Z?-functions which can be de- 
fined operationally as the functions D mmi (R) such 
that, for any rotation R E 50(3), then 

Y^ORn) =^D t mmi {R)Y emi {n) . 

m i 

In this case, it can be shown that the multipole co- 
efficients of the temperature anisotropy map in the 
rotated reference frame are 



D mmi ( R ) a i"i 1 



This expression can be used to compute the correla- 
tion matrix of the aim's in a rotated frame simply as 



(at m a% m ,)R -'^D*J imi (R)D m , m , i (R){ai mi a* l , m , i ) . 

mi 

(A2) 

The Wigner /^-functions take a very simple form 
when we express the rotation matrix R in terms of 
its Euler angles as 

R(a,[3, 1 )=R z (a)-R y {[3)-R z (- / ) . (A3) 

Indeed, for this decomposition we have 

Di m ,(R(a,(3, 7 )) = e« ma+m '^d e mm ,(f3) , (A4) 

where d mm , ((3) = D l mm , (R y (f3)) is a real matrix with 
the following symmetries 

- (-l) m - m 'd m , m (P) , 
dmm'W = d_ m , _ m (p) , 
i-m' ji 



^mm 1 ^ 0) ' ' " i " -111.111 



{-iy- m dt 



dlm'i-P) = (-1) 



m 1 — rn it 



m,m' 



(-1)' 



There exist several explicit and recursive formulae 
to compute these matrices (see A very effi- 

cient recursive procedure can be found in pfij . The 
following formula will be enough to reproduce the 
results presented in this paper. 



18 



2£— 2k-\-m-- m / n \ 2k—m-\-m' 



d l mm >{P) = \[Q- + m)\(l - m)\(l + ml)\{l - m')\ V (-l) fc — \ *' K ,J J , ttt 

mm,\i-> vv /v /v /v ; ^ v i + m - k)\(l - m' - k)\{m' - m + k)\ 



where the sum in k is evaluated whenever the argu- 
ments inside the factorials are non-negative. 

3. Symmetry considerations 

Some consequences of the symmetries of the quo- 
tient manifold on the invariance structure of the cor- 
relation matrix of the ci£ m 's can be deduced directly 
from the transformation rules of the spherical har- 
monic functions under coordinate transformations. 
The results obtained in this way are formal, generic, 
and are very useful in practical computations. We 
end this Appendix by deducing the invariance prop- 
erties the correlation matrix must have, given some 
symmetries of the corresponding quotient manifold. 
These invariance properties have been used in |14| to 
simplify the correlation matrix for the 3-torus, how- 
ever we want to remark here that they are general 
and do not depend on the geometry of the universal 
covering space. 

Let us begin with the invariance properties of 
{aim, a*fj m ') V under rotations around the z-axis. Un- 
der a rotation R z (a) : ip — > ip + a, the function 
Ye m (xi) transforms as 

Y (m (R z (a)n) = e lma Y lm (n) . 

As a consequence the transformation rules for the 
multipole coefficients of a CMB temperature map are 
of the form ag m = e~ lma ae m , and so the correlation 
matrix transforms under this rotation as 

(aim al m ,) T RAa) = e^'~ m ^(a em a% m ,f . (A5) 

We extract two consequences out of IjASjl . First, 
if the quotient space is invariant under a rotation of 
a = 2ir / s around the z-axis, then the correlation 
matrix must be zero unless m = ml mod s. Second, 
if the quotient space is invariant under "any" rota- 
tion around the z-axis, the correlation matrix must 
be zero unless m — ml . In practice, if we take our 
coordinate system such that the fundamental poly- 
hedron of the quotient manifold is oriented so that 
it is invariant under a 2tt/s rotation around the po- 
lar axis, the correlation matrix will present a factor 
^mm' ' an d correspondingly, if the orientation is such 
that the polyhedron is invariant under arbitrary ro- 
tations around the z-axis, the correlation matrix will 
present a factor 5 mm > . 

Let us now take a look at invariance under the 
inversion transformation P : n — * n. Under 
this transformation the spherical harmonic functions 



change as 

Y im (Pn) = (-l) e Y im {n) , 

thus the multipole coefficients ag m change as a£ m = 
(— 1) dim, and as a consequence the transformation 
rule for the correlation matrix is 

(aima* e , m ,)p = (-l) t+i '{ai m al m ,f . (A6) 

Thus, if the fundamental polyhedron is oriented such 
that it appears invariant under the parity transfor- 
mation, the correlation matrix must be zero unless 
£ — £' mod 2, i.e., the correlation matrix will present 
a factor S™ e °, d{2) . 

To end this section, let us consider the reflection 
on the y — plane. This operation changes only the 
azimuthal angle as P y : ip — > —ip, thus the transfor- 
mation rule for the spherical harmonics are 

Y tm (P y n) = Y £ * m (n) , 

the multipole coefficients ag m change as ~ag m — a i m , 
and thus, the transformation rule for the correlation 
matrix is 

(a em a* e , m ,) r Py = (a lm a* t , m ,)* T . (A7) 

It immediately follows that if the fundamental poly- 
hedron is oriented such that it appears invariant un- 
der the reflection on the y = plane, the correlation 
matrix must be real. 



APPENDIX B: CLAUSEN FUNCTIONS 

In this appendix we briefly present some computa- 
tional aspects of the theory of Clausen functions, as 
far as we need them for our purposes. Clausen func- 
tions are periodic functions of period 2tt. There are 
two kinds of Clausen functions, the (^-class and the 
■(/'-class. Clausen (^-functions can be expressed in 
terms of polynomials, while Clausen ^-functions in- 
volve higher transcendental functions, the so-called 
Clausen integrals. Fortunately, we are interested ex- 
clusively in the Clausen ^-functions, thus we will 
develop the details of the theory only for them. The 
Clausen 99-functions are defined as 

00 

<r-^ sin nx 

n—1 

(Bl) 

OO 

, . v-^ cosnx 

n—1 



19 



for s = 1,2,..., and can be calculated recursively 
with the formulae, 



f2a(x) = ((2s) - / tp2s-i(y)dy 



(B2) 



Jo 



where 



is the Riemann Zeta function. 

These recurrence relations are complemented by 
the initial condition 



= E 



sin ni 



= l^2[(2q + l)ir-x] 



(B3) 



e(x-27rg)e(27r(q + l)-a:) 



where 0(x) is the Heaviside step function. Formula 
(|B3|) can be verified by computing the Fourier series 
of the second right hand side. 

Since the Clausen functions are periodic of period 
2tt, we can write 



0{x - 2irq) 0{2ir{q + 1) - x) , 

with fi(x) = 4(tt — x). The recurrence formulae 
(IB2|) yield the following expressions for the Clausen 
functions in the period [0, 2ir], 

f 2s+1 (x) = ^±^—a2( S -r))x 2r+1 + (B4) 



r=0 



(— l) s / irx 



(2s)! (2s + 1)! 



for s = 0, 1, 2, . . . , and 



„2r 



(B5) 



(-1) S / TTX 



„2s-l 



,,2s 



(2a - 1)! (2s)! 



for 8 = 1,2,3,.... 

From the definitions IjBljl we get f2 S +i(n) = 0, 
which can be used to obtain a recurrence formula for 
the Riemann Zeta function of even argument, 



c( 2 *) = E^f^c^--))^- ( B6 ) 



(2s + 1)! 



Writing ((2 s) = g2 S (0)'K 2s , and substituting this into 
i|B6|) we have 



92, 



(°) = E j^+iy. 92{s - r){0) ~ (ferryi 



(-iy s 



The convenience for introducing this notation will be 
apparent in what follows. 

We will now seek for generalizations of the formu- 
lae (IB4|) and (|B5|) . i.e., we look for explicit expres- 
sions for the Clausen functions in the g-th interval 
[27rg, 27r(g + 1)]. Since the Clausen functions sat- 
isfy the periodicity conditions ip2 S -i(2irq) = and 
(p2s(2irq) — ((2s), the recurrence relations (|B2|I can 
be rewritten in the form 



<P2s(x) = C(2s)- / tf2s-i(y)dy 

J2irq 

(B7) 

92s+i(x) = / tp2s{y)dy , 

J2-rrq 

Defining the polynomials f$(x) = f s (x — 2irq), we 
notice that ip s (x) coincides with f§(x) in the interval 
[2irq, 2ir(q + 1)]. This fact, and the expressions 1B7|) . 
allow us to write recurrence formulae analog to (|B2|I 
for the polynomials f§(x) as follows 



fl(x) = 9 2s(qW s - / fl-Mdy, 
Jo 



(B8) 



where 



g2s(q) = 92s{0) + — / fL-Mdy 



(B9) 



I r^q 
92s+l(q) = / f2s(y)dy 



* 2s+1 Jo 



with initial conditions, given by the first Clausen 
function, fl(x) = gi(q)n - § and gi(q) =q+±. 
The expressions (|B8|) can be written in a unified 

way as 



fKx)=9s(q)n s -(-iy / fUMdy- 
Jo 

Using this expression we readily obtain the explicit 
formula, which is the generalization of l|B4|) and (IB5II 
we were looking for, 

/.'(*) = E^i g s -r(q)ir s - r x r - (BIO) 



r=0 



2 U 



20 



where 



where 



H(r, s) 



r l + (-l) s 
2 4 



and L^J is the floor function of x, i.e., the largest 
integer smaller than x. 

The expressions <|B9f) can also be written in a uni- 
fied way as 

/ >27r '? 

ffs ( 9 )= 5s (0) + i-L / /.«_!(!,) dy, 



3,(0) 







if s is even 
if s > 1 is odd 



From this we get the expression analogous to (|B10I 



g s (q) = ffs (0) + 



EC- 1 ) 



""- 1 '- 1) ^s.-r(g)? r -(-i)" ( - M) ^ 



(Bll) 



The polynomials tfo(g) can also be written in the 
canonical form 



(B12) 



where the coefficients are given by A S Q = g s (0), 



A* 



(- i )"]C( -1 )'* (r ~ 1 '"~ 1 i A 



s— r 
n— r 



for < n < 5, and 

"s-l 



r=l 



^ ( _ 1) M(r-l,.-l)l r ^ = r _ 

r_i)A*C*-i.i) ?! 



with initial conditions A J = i and = 1. These 
coefficients are obtained by just introducing i|B12|) 
into (|B11|) and collecting terms. 



APPENDIX C: THE FUNCTION F%,{x) 

In this Appendix we evaluate the function Ffl, (x) 
given by H38|) . We first observe (see Appendix 
that the function Vu' m (x), given by (|4I|I . is an even 
polynomial of (£ + i")-degree. Thus, we begin by 
considering the integral 



1(a) 



P(y) cos ay dy , 



where P(y) is an even analytical function. Integrat- 
ing succesively by parts we get 



1(a) = 2 



sino 



a *■ — ' a 

s=0 



cos a 



(CI) 



E 



ar ' — ' a 

s=Q 



Lg!p(2«+i)(i) 



where P^(x) is the k-th derivative of P(x). 

Making a = nx and P(x) = V Wm (x) in ED , 
substituting I |C1| I in l]38p. and performing the sum 
in n we get 



4£(-l) s 

s=0 



•p(2«) /i \ 
' U'm\ L ) 



r 2s+l 



¥>2 S +i(a:) + (C2) 



V 



(2s+l)/,s 



where (fik(x) is the fc-th Clausen (/^-function defined 
in Appendix iBl 

Since the Clausen functions are periodic functions 
of period 2n, analytic in each period, it follows that 
Ff£, (x) is a piecewise continuous function, analytic in 
each period as well. Thus we will now show how the 
explicit expression for (x) , in the q-th interval 
[2wq, 2ir(q + 1)], given in 1)40(1. comes out. 

Introducing the explicit form for the Clausen ip— 
functions (|B10|) . in the sum of l|C2() yields a huge 
expresion, but a close inspection reveals that it is a 
polynomial in ir/x. The independent term is simply 



e+e' 

2 2- (, + 1)1 



^l(l) 



Vu'm(x) dx 



where the first equality can be deduced by writing 
the Taylor expansion of the integrand of the right 
hand side, and integrating. On the other hand, sum- 
ming up all the coefficients of the (r + l)-th odd 
term, and proceeding as before, we have this term 
equal to 



(-V r ^?i(0)g 2r+1 (q)(l) 



2r+l 
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while the (r + l)-th even term is equal to 

' 7j- ^ 2r+2 

which by the parity of Vu> m {x) is zero. Summing up 
all the terms we finally get iPSJl and 



(-ir« 1) (o) 52 , +2 ( 9 )(^) ; 



APPENDIX D: KNOWN RESULTS FOR 
CLOSED FLAT 3— MANIFOLDS 

In this section we briefly show how we can obtain 
the formulae for the correlation matrix of the aim's 
and the angular power spectrum, currently available 
in the literature, for some closed flat manifolds, as 
well as a simple generalization, i.e., considering the 
observer out of the axis of rotations of the screw 
motions of the covering group. We present explicit 
derivations and formulae for the correlation matrix 
(ai m a\, m ,) , as well as for the power spectrum, in 
order to allow the interested reader to perform their 
own simulations confidently. 

We first give a brief description of flat orientable 
closed 3-manifolds and their covering groups. The 
versions of the diffcomorphic and isometric classifica- 
tions of flat 3-manifolds we present here were given 
by Wolf in [2tJ , and previous descriptions in the con- 
text of cosmic topology were given in |2^ (see 0, 
for alternative descriptions). 

There are six diffeomorphic classes of compact ori- 
entable Euclidean 3-manifolds. The generators for 
the covering groups of the first five classes, Gi — Gs, 
are 71 = (J, a), 72 = (J,b) and 73 = (Ai,c), where 
A\ = I is the identity and 
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for the classes Gi — G5 respectively. It is important 
to remark that these matrices for the rotations are 
written in the basis formed by the set {a, b, c} of 
linearly independent vectors. Thus, the torus Gi is 
generated by three independent translations, while 
for the other manifolds the generators are two in- 
dependent translations and a screw motion along a 
linearly independent direction. The manifold G% is 
the most involved since their generators are all screw 
motions. In the following we present some general 
considerations concerning the classes Gi — G5 only. 

For space forms of the classes G1 — G5, the following 
facts are easily derivable (see [2^ for details): 

1. The vector c is orthogonal to both a and b. 

2. The angle between a and b is a free parameter 
for the class Gi, while its value is fixed to be 
27r/3, it/2 and it/3 for the classes G3, Gi and 
G5 respectively. 




3. Denoting by |c| the length of the vector c, and 
similarly for any other vector, one has that 
|a| = |b| for the classes G3 — G5, while both 
lengths are independent free parameters in the 
class Gi- Moreover, in all classes Gi — G5, |c| is 
an independent free parameter. 

4. Denoting the canonical unitary basis vectors in 
Euclidean space by { e z }, one can always 
write a = |a| e^, b = |b| cos<pe x + |b| sincpe^, 
and c = |c|e z , for the basis {a, b, c}, where <p 
is the angle between a and b, as established in 
the item [2] 

Thus in dealing with manifolds of classes Gi — G5, 
the axis of rotation of the generator screw motion 
can be taken to be the z-axis, and the orthogonal 
part of this generator, in the basis {e x , e y , e z } , is 



(Dl) 



with a — tt, 2tt/3, it/2 and 7r/3 respectively. 
Since the axis of rotation passes through the ori- 
gin, the translational part of the generator 73 is 
c = (0, 0, L z ), where we have put |c| = L z as is 
usual in cosmic topology. 

However, in cosmological applications we need to 
consider the arbitrariness of the position of the ob- 
server inside space. Thus if the axis of rotation is at 
a distance p from the origin (the observer), and its 
intersection with the horizontal plane makes an an- 
gle <fi with the positive x-axis, the translational part 
of the screw motion 73 = (A. c) is 

c = p[cos (j> — cos(0 + a)) e x + 

p[sin (f) — sin(0 + a)] e y + L z e z . (D2) 

In order to perform calculations of the topologi- 
cal signature of CMB temperature maps, we need 
to write the covering group for the manifold under 
study in a compact form. For a torus Gi the prob- 
lem is trivial, since the covering group is generated 
by three independent translations, and thus any two 
isometries commute (see Sec ID II below'), while the 
covering groups for the other closed flat manifolds 
are noncommutative since they contain screw mo- 
tions. 

The generators of the covering groups for the 
classes Gi — G5 satisfy certain relations of the form 



737r i 72 2 



7r i 7 2 " 2 73 , 



where n±, ri2, mi, mi S Z, and they hold whether the 
axis of rotation passes through the origin or not. It 
follows that a generic isometry can always be put in 
the form 



7 



ni no ns 

7i 7 2 7 3 , 



(D3) 



with 7™ 1 = (i",ma), 72 2 = (J,n 2 b), and 
7 3 " 3 =(A h ,n 3 c\\+O h cx) , 
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where A is given by I|D1|) . a = 2ir/s, n 3 = sq + h, 
with q and h integers such that < h < s, the 
parameter s being 2, 3, 4 and 6 corresponding to C? 2 ; 
Qs, Gi and G5 respectively, 



h-l 



°h = J2 Aj 



C|| = L z e zi and c_l = p[cos(f> — cos(0 + a)]e x + 
p[sin 4> — sm(<j) + a)] e y . 

It is now straightforward to compute both the cor- 
relation matrix (o^m&|/ m <) an d the angular power 
spectrum Cg. They all have a simple structure. 
We first describe the general procedure for obtaining 
these expressions and present the results in a unified 
form. We finally specify each case separately. Note 
that, due to (|D2|) . in all of our calculations we are 
considering that the observer may be off an axis of 
rotation of the screw motions of T. 

Upon introducing (|D3|I into 1251) , we transform the 
series of exponentials in a series of Dirac's delta func- 
tions by using l|32[) . The integration of 123f) is then 
immediate in Cartesian coordinates, the general re- 
sult being 



{ag m a e ,. 



(47r) 

V 



Y irm ,(n )Y; m (n )f^(2w0), 

where V is the volume of the manifold, and Z 3 = 
Z 3 \ (0, 0, 0), since the term corresponding to p = 
represents a constant perturbation, and thus is ne- 
glected. 

The function 



1 



L0 —hm e —ik.-O h c 

h=\ 



where oj s is the first complex sth root of unity, is 
a complex modulating term characteristic of the ge- 
ometry and topology of the spatial section of the 
universe model, and depends only on the screw mo- 
tion generators. The vector /3(p) comes from the 
discretization of the wavevector k due to the Dirac's 
deltas (each 2tt(3 is an eigenvalue of the Laplacian 
operator), and is the unit vector in the direction 

of /?. 

Using the property {ai m a% m ,) = {at> m >a* im )* one 
can easily show, by resumming the series, that the 
variances of the multipole moments can be put in 
the general form 



(lair, 



V 



E p 

pGZ 3 



^3 *«(2tt/3) x 



\Ye, 



(n / j)| 2 5R(/r(2 7 ? / 3) 



where 5ft stands for the real part of a complex num- 
ber. The angular power spectrum is then 



= Y E i*«(2^)S,(2^) 



where 



EAk) =2lTT Z l^ m K)| 2 5ft(.^(k)) 



m= — I 



can be evaluated using the Addition Theorem for 
Spherical Harmonics yielding 



5*(k) = - 

s 



s-l 



1 + Pi(cos8 Kh ) cos(k • O h c) 



h=l 



where 



cos = cos 9k + sin 6*k cos • 



2rh 



1. Rectangular torus Q\ 

The generators for the rectangular torus are the 
translations a = L x e Xl b = L y e y and c = L z e z , 
thus a generic isometry of its covering group can be 
written as 7 = (I, r), with r = n x L x e x + n y L y e y + 
n z L z e z , and n x ,n y ,n z € Z, i.e., the covering group 
of Q\ is parametrized by Z 3 . 

It follows immediately that, for a rectangular 
torus, the expression (|25|) takes the form 

n(EZ 3 

= (2^) 3 ^ 6(k x L x -2ir Px ) x 
pez 3 

8{kyL y -2-Kpy) 8{h z L z 2 , Kp z )Yi m (nit) . 

Following the procedure described above one gets 
/£(k) = 1, and (3 X = and 0, = £. In 

particular, we have the well known result 

c < = y E i*«( 2 ^)- 

P ez 3 

2. Rectangular t/2 

The generators for the rectangular Q2 are 71 = 
(I, a), 72 = {I, b) and 73 = (A,c), with a = L x e x , 
h = L y e yi c = 2pcosc/)e x + 2psm(j)e y +L z e z , and A 
given in 1|D1|) with a — ir. They satisfy the relations 
7i737i = 73 and 727372 = 73, which allow to write 
any isometry of the covering group by l|D3(l with 



„ 3 _ f (J,n 3 C||) if n 3 is even . . 

73 ~ \ (A,n 3 cn +c ± ) if n 3 is odd ' [ ' 



where cj_ = 2pcos<pe x + 2psin</>e a . 

It follows from (|D3|) and 1D4(1 that the expression 
(|25|l takes the form 



pS2 
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T L( k ) = e- i{ - n * k * L * +n -» k v L * +2n * h » L *"> [l + (-l) m e" lk - c ] y /m (n k ) 

= (2?r) 3 Y s (k x L x - 2np x )S(k y Ly - 2irp y ) S(k z L z - np z ) Yg m (n k ) /£(k) , 



where we have put n\ = n x , ri2 = n y , and 71-3 = 2n z 
or 2n z + 1, depending on whether 713 is even or odd. 
The components of f3 are /3 X = jf-, (3 y = jf-, and 

^* ~ 2L 2 ■ 

3. g 3 

The generators for a manifold of class Q3 are 71 = 
(J, a) j 72 = (I,b), and 73 = (A,c), with a = Le^, 
b = — -j (e x — v^e^), c = I (3cos0 + v^sin^) e x + 
I (3 sin — y3 cos (/>) e y + L z e z , and A given in (|D1|) 



I 

with a = 27r/3. They satisfy the relations 7 2 ~ 1 737i = 
73 and 71727372 = 73, which allow us to write any 
isometry of the covering group by l|D3|) with 

!(/, ri3C||) if 113 = mod 3 

(A, ri3C|| + cj_) if 713 = 1 mod 3 , 
(A 2 , n 3 C|| + 09C ± ) if n 3 = 2 mod 3 

(D5) 

where cj_ = § (3cos</> + V3 sin <j>) e x + | (3sin</> — 
V^cos 0) e H . 

It follows from (|D3|) and 1D5(1 that the expression 
l|25[) takes the form 



_ ^""^ iJrixfei£+™ B (T^' c B — ifecj L+3n z k z L z ^ j- 



1 — m — ik-c 1 —2m — ik-Oocl /„ \ 

1+uj 3 e +uj 3 e I Y tm (n k ) 



2tt) 3 ^ 5(fc x L-2 7 rp :c )j([^fc J/ -iA: x ]i-27rp J/ ) 5 (k z L z - %p z ) Y im (n k ) fa (k) , 



and n 3 = 3n 2 



where we have put m — n x , n-i = 
3n z + 1 or 3n z + 2 according to |D5j). We also get 

(3 X = %,0y = $ (2py+p x ), and B z = 



The generators for a manifold of class G4 are 
7i = (J) a), 72 = (I,b), and 73 = (A,c), with 



a = Le x , b 
p (sin 4> — cos (/i 



Le,, 



p (cos (/> 



e y + L z e z , and A given in I|D1|) with 



a = 7r/2. They satisfy the relations 7 2 1 



737i 



73 



and 717372 = 73, which allow us to write any isom- 
etry of the covering group by (|D3|I with 



1^3 
7 3 



(7, rt 3 cy) 
(A, 713C11 - 
(A 2 ,n 3 cy 
(A 3 ,n 3 C|, 



■c±) 
f 2 c x ) 
f 3 c±) 



if 713 = mod 4 
if 713 = 1 mod 4 
if «3 = 2 mod 4 
if «3 = 3 mod 4 



(D6) 

where cj_ = p (cos + sin <j)) e x + p (sin ip — cos (/)) e y . 

Similarly using (|D3|) and ljD6|) , the expression (|2*5)l 
takes the form 



T L( k ) 



-i{n x k x L+n y k y L+An z k z L z ) 



-t , \ ^ —hm - 

1 + > w 4 e 



(2^) 3 J] 5(fc x L - 2^ x ) S(k y L - 2tt Pv ) 6 (k z L z - § p a ) F, m (n k ) /£(k) , 



r 



where we have put n\ = n x , n-i = n y , and 713 = 4n 2 , is the first complex 4th-rooth of unity. We also 

4n z + 1, 4n z + 2 or 4n z + 3 according to ((D6|) . and get /3 X = 2^, (3 y = and /3 Z = 
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The generators for the rectangular Q§ are 71 = 
(J, a), 72 = (I,b), and 73 = {A, c), with a = Le x , 
b = ■§■ (e x + v^e^), c = § (cost/) + V3sin</>) e x + 
£ (sin^ — \/3cos0)e y + L 2 e 2 , and A given in (|D1|) 
with a — 7r/3. They satisfy the relations 7 2 ~ 1 737i = 
73 and 7i7 2 _1 7372 = 73, which allow us to write any 
isometry of the covering group by l|D3(l with 



where cj_ = § (cos 4> + V3 
\/3cos 0) e„. 



sine 



(sin < 



Tlx 

7 3 



(I,«3C||) 

(A,n 3 C|| - 
(A 2 ,n 3 C|| 

(A 3 ,n 3 C|| 

(A 4 ,n 3 C|| 
(A 5 ,n 3 C|i 



•c±) 
f0 2 ci) 
f 3 c±) 
f0 4 ci) 

f 0BC±) 



if 713 = mod 6 
if n 3 = 1 mod 6 
if rj-3 = 2 mod 6 
if tt-3 = 3 mod 6 
if n 3 = 4 mod 6 
if 713 = 5 mod 6 



(D7) 



Using (|D3|> and i|D7|) , the expression (|25(l takes the 
form 



rL(k) = ]T 



1 1 \ ^ —hm 

1 + 2^ UJ 6 e 



Y« m (n k ) 



h=l 



(2tt) 3 6 ( k * L - 2 *Px) 6 ([^ k v + \ k *\ L ~ 2 ^Pv) 6 { k * L * ~ f Vz) ^m(n k ) /„(k) , 



where we have put m — n x , n 2 — n y , and = 6n z , get [3 X = (3 y = ^ (2p y — p x ), and (3 Z — J^-. 
6n z + 1, 6n z + 2, . . . , 6n 2 + 5 according to (|D7jl . and 
W6 is the first complex 6th-rooth of unity. We also 
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